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An unusual quasi-two-dimensional heavy band mass van Hove singularity (vHs) lies very near the 
Fermi energy in MgCNis, recently reported to superconduct at 8.5 K. This compound is strongly 
exchange enhanced and unstable to ferromagnetism upon hole doping with ~12% Mg^Na or Li. The 
1/4-depleted fee (frustrated) Ni sublattice and lack of Fermi surface nesting argues against competing 
antiferromagnetic and charge density wave instabilities. We identify an essentially infinite mass along 
the M-F line, leading to quasi-two-dimensionality of this vHs may promote unconventional p-wave 
pairing that could coexist with superconductivity. 



The discovery of ~40 K superconductivity in MgB2 0] 
has spurred interest in searching for superconductivity in 
unUkely materials, and other discoveries are uncovering 
previously unanticipated relationships between FM and 
superconducting states. Of this latter category, there are 
now several examples, such as the magnetic organometal- 
lic (BETS)2FeCl4, where superconductivity is actually 
induced Q by a strong applied magnetic field rather than 
being destroyed by it, and the intermetallic UGe2, where 
superconductivity occurs ||^,^ in spite of strong ferromag- 
netism (FM) and coexists with it to the lowest tempera- 
tures, and is almost certainly triplet paired. The question 
of FM on the surface of the exotic superconducting 
oxide Sr2Ru04 is yet another aspect of the strong rela- 
tionship between FM and superconductivity that is yet 
to be understood. 

A new compound that would not be considered a likely 
candidate for superconductivity is MgCNis, whose con- 
duction electrons are primarily derived from Ni which is 
itself a FM, yet it superconducts at 8.5 K. ^ This per- 
ovskite compound (a = 3.812 A) can be regarded as fee 
Ni that is expanded by 8%, one quarter of the Ni replaced 
by Mg, then C atoms put into the octahedral sites. Par- 
tial replacement of Ni by Co or Cu reduces (or destroys) 
Tc, 1^ while other substitutions have not been reported. 
The Hall coefficient is hole-like but strongly dependent 
on temperature, unlike a conventional metal (Fermi liq- 
uid), while the measured critical field Hc2 has a conven- 
tional shape. ^ Tunneling spectra show a strong zero 
bias anomaly that has been interpreted as evidence of 
strong coupling superconductivity of an unconventional 

type. 

In this paper it is shown that this compound, in addi- 
tion to being superconducting, is also an incipient FM, 
which can be driven to ferromagnetism by partial («12%) 
replacement of Mg with a monovalent metal such as 
Li or Na. This close proximity of superconductivity to 
magnetism itself implies unconventional pairing, and our 
identification of a van Hove singularity (vHs) with quasi- 
two-dimensional (2D) character provides further support 



for unconventional pairing and possible coexistence of 
FM and superconductivity, two types of collective order 
that are usually antagonistic. 

The perovskite structure itself is rather unusual for 
such an intermetallic compound, since perovskites much 
more commonly have a strongly negative ion (0^~ or a 
negatively charged halide) on the site occupied by Ni in 
this compound. It is essential first to understand the 
character of the charge carriers, for which purpose we 
have carried out full potential, all-electron density func- 
tional based calculations. The resulting spectral dis- 
tribution of the electronic states (for the experimental 
lattice constant a=3.812 A) is shown in Fig. 1, and is in- 
distinguishable from that presented by Singh and Mazin. 
Ip^ The conducting states at the Fermi level are dom- 
inated by Ni dxz,dyz and also d^-2_y-2 character, in the 
local coordinate system in which the z axis is directed 
toward the two neighboring C atoms. The remarkable 
feature of this compound is the sharp peak in the den- 
sity of states (DOS) just 45 meV below the Fermi level. 
p3[ This peak results from a van Hove singularity (vHs) 
arising from a remarkably fiat, primarily Ni 3d derived, 
band at and around the M = (l,l,0)7r/a point in the 
simple cubic Brillouin zone. Thus it is high mass Ni 
3d holes that form the superconducting pairs, and Ni-Ni 
hopping is important for their transport. 

To understand more clearly the origin of this peak, a 
fictitious material D^+CNis was studied, i.e. Mg was re- 
moved but its two valence electrons were retained. While 
the C atom does have an appreciable effect on the bond- 
ing and the resulting density of states (not pictured) , for 
practical purposes the Mg simply gives up its to valence 
electrons to the bands (formed mainly by Ni) and has al- 
most no other effect, as can be seen in Fig. 1. We utilize 
this important point below. If there were 0.5 electron 
less per cell, Ej? would lie just at the peak in the DOS 
where N(E) is a factor of eight larger. 

In addition to promoting superconductivity, a large 
value of N(Ef) (equal to 2.4 states/eV-spin here) leads 
an exchange-enhanced magnetic susceptibility x that 
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strongly opposes singlet superconductivity, or possibly 
even to a FM instability (where % — > oo) which is in- 
compatible with singlet superconductivity. This latter 
scenario applies to Sr2Ru04, [Q which is a nearly FM 
superconductor (but only below 1.5 K) and is now un- 
derstood to be a parallel-spin-paired (triplet) supercon- 
ductor. Density functional calculations are very reliable 
in calculating this tendency toward magnetism, and in- 
deed the instability to FM, particularly in intermetallic 
compounds such as MgCNia. The enhanced (observed) 
susceptibility p5|] is given by 



X = 



Xo 



1 - N{Ef)I 



SXo, 
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where Xo — 2/i^N(Ei?) is the bare susceptibility obtained 
directly from the band structure and / is the exchange 
interaction. 

We have calculated / « 0.29 ± 0.01 eV in two ways. 
One, which demonstrates directly our main thesis that 
MgCNia is close to ferromagnetism, was a calculation for 
ordered Mgi/2Lii/2CNi3, which effectly simply removes 
0.5 valence electron from the cell. This material is pre- 
dicted to be FM, and the exchange splitting A^x between 
majority and minority bands (Fig. 2) gives /= 0.30 eV 
from the relation Ae^; = Im where m is the FM mo- 
ment in units of /i^. The other calculation of / resulted 
from fixed spin moment calculations, in which the 
energy E{m) is calculated subject to the moment be- 
ing constrained to be m. The behavior at small m is 
E{m) = (l/2)x~^»Ti^ from which / — 0.28 cV can be 
extracted from Eq. (1). Singh and Mazin |l^ have ob- 
tained a similar value for the Stoner parameter. Thus 
5=3.3, and it is certainly unexpected for a conventional 
(singlet) superconducting state to survive so near a FM 
instability, especially when the superconducting carriers 
are heavy and are the same ones that will become mag- 
netic. 

To quantify how near this system is to being FM, we 
have carried out (i) a series of virtual crystal calcula- 
tions for Mgi_i-NaA-CNi3 (justified by the results shown 
in Fig. 1) to find the concentration 5cr of the FM critical 
point, and (ii) an extended Stoner analysis. The two 
results are consistent in predicting the onset of FM at 
dcr — 0.12. The ordered magnetic moment m{S) versus 
hole doping level is shown in the inset in Fig. 3, where 
it is evident that, in the absence of superconductivity, a 
moment grows as m{S) — Q{S — ScrY^'^ for small S — 6cr 
beyond the critical concentration. The behavior of m{6) 
in the small m limit can be obtained analytically from 
an expansion of the DOS N{m) averaged over the states 
within ±m of Ep: 
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and using iV(0, Scr) = 1// to obtain the square root law, 
with 
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There are experimental indications from tunneling [Tc| ] 
that the superconductivity in MgCNia may arise from 
triplet pairing. |l^ Our results show that hole doping 
with Na or Li will be an excellent way to probe this pos- 
sibility. While it is unexpected that singlet supercon- 
ductivity would occur at all in a Ni compound that is as 
strongly exchange enhanced as MgCNia is (the increasing 
enhancement that diverges as (S — > Scr should kill singlet 
superconductivity very quickly) , for triplet pairing the in- 
creasing magnetic correlations may provide the coupling 
(as recently argued for heavy fermion superconductors 
|19| ) and enhance Tc as (5 — > Scr- 

We now consider in more detail the electronic struc- 
ture and its implications, especially in the critical region 
S « 0.12. In this region is still 35 meV from the 
vHs, so (near) divergence of N(e) per se is not a crucial 
consideration. However, FM, AFM, and pairing suscepti- 
bilities should be considered, and they will involve energy 
denominators Se = ^ — £j_|_q-„, where n,TO label the 
three M points where the vHs occur. For n = m, Se will 
be small for small q, which relates to the FM instability 
we have considered above. Since there are vHs at three 
inequivalent M points, there will also be a potential AFM 
instability near a vector Q that spans two inequivalent 
vHs. These values of Q are in fact equal to M, so AFM 
tendencies are peaked at Q = (1, l,0)7r/a and symmetri- 
cally related points. 

The dominant instability is determined largely by 
phase space availability. The band giving the vHs is quite 
flat (to within 50 meV) in a roughly cubic region of side 
Tr/a centered on each M point, which totals to 3/8 of the 
Brillouin zone volume. Band masses at the vHs are quite 
large: m* « 15 along M ~ X, m* k -15 along M - R. 
Along M — r, however, this band is constant to within 
0.2 meV for a distance of « i-^, i.e. it has effectively 
an infinite effective mass along this line. The dispersion 
relation in the X — M — V plane is pictured in Fig. 4. 
Since the plane containing this dispersionless line is ori- 
ented differently for each M point, there is no pronounced 
nesting of the vHs. However, because of the lack of dis- 
persion along one direction, the shape of the N(E) peak 
in Fig. 1 has the logarithmic divergence characteristic 
of 2D models. This vHs is remarkably insensitive to vol- 
ume, being essentially unchanged relative to Ej? when the 
lattice constant is decreased by 5%, which makes pres- 
sure a relatively ineffective tool (except perhaps to tune 
the Curie temperature beyond the critical point). 

The structure of the vHs suggests that the resulting 
physics could be a three dimensional generalization of 
the 2D vHs models. The strong exchange enhancement 
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and the proximity to a FM instability is most consis- 
tent with triplet pairing, appropriate forms of which can 
coexist with FM as in UGe2, or with incipient FM as 
in Sr2Ru04. There are numerous possibilities of p (or 
higher odd) symmetry for a cubic system [ pO[ , some of 
which have gaps but many of which have nodes and there- 
fore are gapless. It typically requires high quality ther- 
modynamic and spectroscopic data, which are not yet 
available, to determine the character of the order param- 
eter. 

The theoretical case for superconductivity arising from 
longitudinal magnetic fluctuations was laid out by Fay 
and Appel, who found that Tc should peak near 
^ = ^cr (on either side) but would vanish at the critical 
point. Machida and Ohmi emphasize that a non- 
unitary triplet state is most likely in a case that may ap- 
ply to MgCNis, since it intrinsically breaks time- reversal 
symmetry as does ferromagnetism. In such non-unitary 
phases, a magnetic field may enhance Tc. psf 

The most intriguing possibility is that (triplet) super- 
conductivity might coexist with ferromagnetism, as re- 
ported recently for UGea. Mgi_5(Na,Li)5CNi3 
provides the new possibility, if indeed coexistence oc- 
curs due to triplet pairing, of studying the emergence 
of itinerant ferromagnetism (as 5 crosses dcr) within an 
existing superconducting phase. The phenomenological 
theory of coexistence in just such a case has been put 
forward recently, concluding that the heat capacity 
has a linear-in-T term that is strongly dependent on the 
magnetization. Solutions of the Eliashberg equations for 
a spin fluctuation system near the quantum critical point 
p3| suggest that triplet superconductivity might not be 
as strongly favored near the critical point as might have 
been anticipated. Hole-doped MgCNia appears to be an 
excellent system to use as a probe of these fundamental 
questions. 

The FM instability a,t S = Scr was demonstrated above. 
Can AFM ordering of the Ni atoms, such as is encour- 
aged by inequivalent vHs in other systems, occur here? 
The Ni sublattice is a 1/4 depleted fee lattice, but the 
periodic depletion does not alter the well known frustra- 
tion of Neel order on the fee lattice. Thus while short 
range AFM correlations may be strong, AFM ordering 
will be opposed by frustration. 

To summarize, we have shown that the superconductor 
MgCNis is near a ferromagnetic instability that can be 
reached by hole doping on the Mg site. The effective car- 
riers are Ni-derived holes of very high band mass (likely 
enhanced by dynamic spin fluctuations and phonons). 
The FM instability is related to an unusual quasi-2D 
heavy mass van Hove singularity less that 50 meV be- 
low E^. This quasi-2D character supports earlier sug- 
gestions that its superconductivity is unconventional in 
nature, and we suggest that hole-doping is an ideal way 
to probe the onset of ferromagnetism in the supercon- 
ducting state. 
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FIG. 1. The density of states of MgCNis (top panel), 
showing the nearly filled Ni d states and the extremely sharp 
and narrow peak just below Ef arising from a van Hove sin- 
gularity. The inset gives a blowup of the peak, with the 
shaded portion indicating the dominant Ni 3d contribution. 
The lower panel shows the density of states for D^^CNis (see 
text) . 



12 




-3.5 -2.5 -1.5 -0.5 0.5 1.5 

Energy (eV) 



FIG. 2. Ferromagnetic states 

obtained from self-consistent calculation. Density of 
states of ferromagnetic S = 0.5 materials. Top panel: the 
ordered compound Mgo.sLio.sCNis. Bottom panel: a virtual 
crystal result for Mgo.sNao.sCNia. This level of doping results 
in a filled majority van Hove peak and an empty minority van 
Hove peak. The magnetic moments are 0.83 /ib and 0.95 /is, 
respectively. 
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FIG. 3. The mean value N(m) of the density of states 
around Ep necessary to create a magnetic moment m, ver- 
sus m. Hole-doping concentrations < S < 0.45 are shown. 
The dashed line indicates l/I (1=0.29 eV). Solid curves give 
results from a rigid band treatment based on the MgCNis 
{5 = 0) DOS. The inset gives the predicted value of the fer- 
romagnetic moment versus the hole concentration from the 
Stoner model (solid line) and from specific self-consistent vir- 
tual crystal calculations (dots) for Mgi-jNa^CNis, which in- 
dicates the consistency. The square gives the moment for the 
ordered compound Mgi/2Lii/2CNi3 discussed in the text and 
in Fig. 2. 
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FIG. 4. Surface plot (and contour plot below) of the van 
Hove singularity in £fc (relative to Ef) in the T — M — X 
plane, with M at the right hand corner (planar coordinates 
are given in units of 2-k ja). Note the extreme flatness of the 
dispersion along the diagonal M ^ F. The negative effective 
mass sheet of the vHs lies in the third direction. 
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